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The quark-gluon plasma (QGP) generated in relativistic heavy-ion collisions could be locally parity-odd. In 
parity-odd QGP, the electric field may induce a chiral current which is called the chiral electric separation effect 
(CESE). We propose two possible observables for CESE in Cu -l- Au collisions: The first one is the correlation 
fag = <cos[2(0„ + 4>g — 2'Prp)]); the second one is the charge-dependent event-plane angle “Pj with q = ± being 
charge. Nonzero Af = - f^ame and A*? = <|'PJ - T'jl) may signal the CESE in Cu -l- Au collisions. Within a 

multiphase transport model, we study how the final state interaction affects these observables. We find that the 
correlation jag = <cos((^„ + fg — 'Prp)) is sensitive to the out-of-plane charge separation caused by the chiral 
magnetic effect and to the in-plane charge separation caused by the in-plane electric field, but it is not sensitive 
to the CESE. On the other hand, Af and AT* are sensitive to the CESE. Therefore, we suggest that the future 
experiments measure the above observables in Cu-tAu collisions in order to disentangle different chiral and 
charge separation mechanisms. 

PACS numbers: 25.75.-q, 12.38.Mh, 11.30.Er 


I. INTRODUCTION 

Relativistic heavy-ion collisions generate not only ex¬ 
tremely hot quark-gluon plasma (QGP) but also extremely 
large magnetic fields due to the fast motion of the colliding 
ions. Recent detailed calculations revealed that the maximum 
magnetic fields in Au + Au collisions at energies currently 
available at the BNL Relativistic Heavy Ion Collider (RHIC) 
can reach Sm^ ~ 10'* Gauss while in Pb + Pb collisions at 
energies available at the CERN Large Hadron Collider (LHC) 
they can reach 60m^ ~ 10*® Gauss [1-5]. Under such large 
magnetic fields, some novel quantum phenomena can possi¬ 
bly happen. The most intriguing ones are the so-called chi¬ 
ral magnetic effect (CME) [6-10] and chiral separation ef¬ 
fect (CSE) [11, 12]: They occur in parity- and charge-odd 
regions in QGP and can result in charge and chirality sepa¬ 
rations along the direction of the magnetic field, respectively. 
Recent experimental measurements of the charge azimuthal 
correlation [13], 

yag = <cos(0„ + 2'1'rp)), (1) 

where tpa and (pg are the emission azimuthal angles of parti¬ 
cles with charges a and (3 and 'Prp is the reaction plane angle, 
showed some consistent features with the expectation of CME 
at both RHIC and LHC energies [14-16]. Another important 
experimental test of CME and CSE is the observation of the 
charge asymmetry in the elliptic flow of pions [17, 18] that 
is consistent with the expectation of a chiral magnetic wave 
(CMW) [19, 20] — a collective mode arising due to the inter¬ 
play between CME and CSE in the presence of the magnetic 
field. However, one must notice that there are other back¬ 
grounds that contribute to the experimental observables, see 
Refs. [21-29] for discussions. 
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Heavy-ion collisions can also generate strong electric fields 
due to the event-by-event fluctuation of the proton positions 
in the ions [3, 4] or due to asymmetric colliding geometry 
(for example, in Cu + Au collisions [30-32]). It was pro¬ 
posed that the electric fields may also induce chiral current 
and chiral separation in QGP, which is called the chiral elec¬ 
tric separation effect (CESE) [33-36]. It was also proposed 
that the Cuh-Au collisions may provide us a good chance to 
detect CESE because there is a strong electric field directing 
from Au to Cu due to the charge asymmetry between Au and 
Cu nuclei [33]. Recently, the electromagnetic fields in Cu H- 
Au collisions were studied in details, and it was found that the 
large in-plane charge separation effect (in-plane CSE) induced 
by the strong in-plane electric fields could strongly suppress 
or even reverse signs of jag [30]. 

In this paper, we propose two observables that are designed 
solely for the detection of CESE. In addition, it is crucial to 
take the final state interactions into account for any model cal¬ 
culations in order to link the initial anomalous transports to the 
experimental data since heavy-ion collisions undergo compli¬ 
cated dynamical evolutions which involve many final interac¬ 
tions. We study the effects of the final state interactions by 
using a multi-phase transport (AMPT) model which success¬ 
fully describes the main evolution stages of heavy-ion colli¬ 
sions. By introducing appropriate initial dipolar or quadrupole 
charge distributions, the AMPT model can successfully de¬ 
scribe both the charge azimuthal correlation jag [37, 38] and 
the charge asymmetry of the pion elliptic flow [39] in Au + Au 
collisions at the top RHIC energy. In this work, we introduce 
different kinds of initial charge separations which are used to 
mimic different initial chiral effects into the initial condition 
of the AMPT model, and predict some observables which can 
be used to test whether these effects can be observed in the 
final state of Cuh-Au collisions at = 200 GeV. 

The paper is organized as follows. In Sec. II we setup the 
numerical simulations, in Sec. Ill we show our main results. 
Einally we summarize and discuss in Sec. IV. Throughout this 
paper, we use the natural units Ti - kg - c - \. 
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II. GENERAL SETUP 
A. AMPT model 

The AMPT model with a string melting scenario is utilized 
in this work [40]. The AMPT model is a dynamical trans¬ 
port model which consists of four main components: initial 
condition, parton cascade, hadronization, and hadronic rescat¬ 
terings. The initial condition, which includes the spatial and 
momentum distribution of participant matter, minijet partons 
production and soft string excitations, is obtained through the 
HIJING model [41,42]. The parton cascade starts the partonic 
evolution with a quark-anti-quark plasma from the melting of 
strings. Parton scatterings are modelled by Zhang’s parton 
cascade (ZPC), which currently only includes two-body elas¬ 
tic parton scatterings using cross sections from pQCD with 
screening masses [43]. The parton interaction cross section 
is set as 10 mb in our work, by which the model has shown 
good abilities to describe many key experimental observables 
at RHIC [44-49]. A quark coalescence model is then used 
to combine partons into hadrons when the system freezes out. 
The evolution dynamics of the hadronic matter is described by 
a relativistic transport (ART) model [50]. Because the current 
implementation of the ART model does not conserve the elec¬ 
tric charge, we in this study consider only resonance decays 
and hadronic scatterings are switched off to ensure the charge 
conservation. 


B. Introducing the initial charge separations into the AMPT 
model 

Let us consider the chiral magnetic effect (CME) in Cu 
+ Cu collisions first. This can serve as a test of the AMPT 
model. The coordinate system is set up so that the x-axis is in 
the reaction plane, i.e., T'rp = 0, and the y-axis is perpendic¬ 
ular to the reaction plane with the z-axis being the direction 
of the projectile nucleus. To mimic the CME, we introduce 
an initial charge separation into the initial state of the AMPT 
model, since the charges are not separated but distributed ran¬ 
domly in the normal AMPT model. To separate a percent¬ 
age of the initial charges, we follow the procedure of a global 
charge separation scenario, which has been employed in Ma 
and Zhang’s previous work [37]. We switch a percentage of 
the downward moving u quarks with those of the upward mov¬ 
ing u quarks in such a way that the total momentum is con¬ 
served, and likewise for d and d quarks, where the percentage 
should presumably depend on the impact parameter b or cen¬ 
trality, because the magnitude of the averaged magnetic field 
is b dependent. 

In Pig. 1 we show our AMPT results of the charge az¬ 
imuthal correlation jap (from the normal AMPT model and 
the AMPT model with the initial CME-like charge separa¬ 
tion) as well as the experimental data for CuH-Cu collisions 
at = 200 GeV. We find that, to fit the CuH-Cu experi¬ 
mental data, the percentage of initial CME-like charge sepa¬ 
ration /% should be proportional to the impact parameter b 
with a slope of 1.56, i.e. / = 1.56l7/fm. This is consistent 



FIG. 1: (Color online) Centrality dependence of (cos(0j, -l- (j)^)) in 
Cu-tCu collisions at y5[^=200 GeV. The open symbols represent 
the results from the normal AMPT model without anomalous effects 
and the AMPT model with an initial CME-like charge separation, 
and the solid symbols represent the experimental data [14]. 


with the fact that the averaged magnetic field is proportional 
to b [3, 4]. Comparing with the normal AMPT result, the 
AMPT result with the initial CME-like charge separation can 
well describe CuH-Cu data measured by the solenoidal tracker 
at RHIC (STAR) experiment. 

We now turn to study the Cu + Au systems. The co¬ 
ordinate system of Cu + Au is set up similarly to that of 
CuH-Cu, but the Au nucleus is set to the left of the Cu nu¬ 
cleus, i.e., the direction of the total electric field is rightward. 
Eigures 2 (a)-(f) show the net electric charge distributions in 
the transverse momentum space of the initial partonic states 
in the AMPT models with different kinds of initial charge 
separations for centrality bin of 30-40% in Cuh-Au collisions 
ys^^=200 GeV. After introducing the CME effect into the 
normal AMPT model, the net electric charge distribution is 
changed from that in Pig. 2 (a) to that in Pig. 2 (b). Due 
to the strong electric fields in the in-plane direction in Cu H- 
Au collisions, there may appear also the in-plane charge sep¬ 
aration effect (in-plane CSE) [30]. To simulate the in-plane 
CSE, we take an analogous setup as above by switching the 
values of a percentage of the leftward moving u quarks 
with those of the rightward moving u quarks, and likewise 
for d and d quarks. The b dependence of the percentage is 
assumed to be the same as that for the out-of-plane separa¬ 
tion. The corresponding net electric charge distribution is 
shown in Pig. 2 (c). Now we take the chiral electric sepa¬ 
ration effect (CESE) into account. We consider the situation 
where the CESE, CME and in-plane CSE happen simultane¬ 
ously (denoted as CESE-nCME-Hin-plane CSE) [33]: Por the 
quarks with py > Owe switch a percentage /% of the leftward 
moving u quarks with those of the rightward moving u quarks 
(likewise for d and d quarks); while for quarks with p,. < 0 
we switch half of the percentage (i.e. 0.5/%) for the leftward 
moving u (d) quarks with those of the rightward moving u (d) 
quarks. Por the initial charge separation /%, we apply the 
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FIG. 2: (Color online) The net electric charge distributions in the transverse momentum plane in the initial partonic states for different settings 
of AMPT models for centrality bin of 30-40% in Cu-tAu collisions yj^^=200 GeV. (a) Normal setting without any anomalous effect, (b) 
AMPT with out-of-plane charge separation which mimic the chiral magnetic effect (CME). (c) AMPT with in-plane charge separation effect 
(in-plane CSE) caused by the in-plane electric field, (d) AMPT with in-plane CSE plus a quadrupolar distribution due to chiral electric 
separation effect (CESE) and CME. (e) The same as (d) but with the strength of CME doubled, (f) The same as (d) but with the strength of the 
in-plane CSE doubled. 


same /t-dependent initial charge separation as that of Cu-i-Cu 
collisions at yi^^=200 GeV. In this case, this It-dependent 
percentage is expected to provide a lower limit for the initial 
charge separation percentage, since the Cu-i-Au system has a 
little bit stronger magnetic field than the Cu-i-Cu system at the 
same impact parameter [30]. The net effect is equivalent to 
a configuration with a in-plane dipole plus a quadrupole as 
illustrated in Fig. 3. The net electric charge distribution for 
CESE-i-CME-i-in-plane CSE is shown in Eig. 2 (d). To study 
the case where the CME is much stronger, we do the CME 
switching once more after the above CESE-i-CME-i-in-plane 
CSE switching, which is shown in Eig. 2 (e). If the in-plane 
CSE effect is much larger, we do the in-plane CSE switch¬ 
ing once more after the CESE-i-CME-i-in-plane CSE switching, 
see Eig. 2 (f). In these ways, we have six different kinds of 
initial charge distributions which mimic different initial chi¬ 
ral effects. We use them as the different inputs for the initial 
condition of the AMPT model, and then we extract various 
hadronic observables to test whether the initially embedded 
effects can survive after the final state interactions. The re¬ 
sults are presented in the following. 



FIG. 3; (Color online) Illustration of the charge configuration of 
chiral magnetic effect (CME) plus chiral electric separation effect 
(CESE) plus in-plane charge separation effect (in-plane CSE) in Cu 
+ Au collisions. 


III. RESULTS AND DISCUSSIONS 
A. Correlation 

We first consider the correlation jap as defined in Eq. (1). 
The centrality dependence of jap in Cu-i-Au collisions at 
ys^^=200 GeV from different settings of the AMPT model 
are shown in Pigs. 4 (a)-(f), where the centrality bins are de¬ 
fined by using different ranges of impact parameter. Also 
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FIG. 4: (Color online) Centrality dependence of = (cos(i^a. + 0 ^ 5 )) from different initial settings of AMPT models (open symbols) in Cu+Au 
collisions at yj^^=200 GeV. Also shown are the experimental data for Cu+Cu collisions at yi'^^=200 GeV (solid symbols) [14]. Different 
panels are in one-to-one correspondence with Fig. 2. 


shown is the experimental data from Cu + Cu collisions which 
are used as the baseline for comparison. Figure 4 (a) shows 
that the normal AMPT model without any initial charge sep¬ 
aration gives a zero opposite-charge correlation and a nega¬ 
tive same-charge correlation which have much smaller mag¬ 
nitudes than the experimental data for CuH-Cu collisions at 
yi^^=200 GeV. After introducing the chiral magnetic effect 
(CME) into the AMPT model, as shown in Fig. 4 (b), the 
opposite-charge correlation becomes positive and the same- 
charge correlation becomes more negative, which present sim¬ 
ilar magnitudes as the CuH-Cu data. It is consistent with the 
previous AMPT works about Auh-Au collisions [37] and the 
above results about CuH-Cu collisions. If the initial charges are 
separated along the reaction plane direction [mimicking the 
in-plane charge separation effect (in-plane CSE) caused by 
the strong in-plane electric field], comparing with the CME 
case, opposite-charge and same-charge correlations reverse 
their signs, i.e. a negative opposite-charge correlation and a 
positive same-charge correlation are observed and their mag¬ 
nitudes are comparable to the CuH-Cu data, as shown in Eig. 4 
(c). This result is consistent with recent study in Ref. [30] 
and more discussions can be found there. It is also consis¬ 
tent with the charge asymmetry of direct flow vi suggested 
in Ref. [31]. If the chiral electric separation effect (CESE), 
CME, and in-plane CSE happen together in the initial stage 
[Eig. 4 (d)], both the opposite-charge and same-charge corre¬ 


lations are strongly suppressed to the levels close to the nor¬ 
mal AMPT model [Pig. 4 (a)], which indicates the three ef¬ 
fects almost cancel out. This can be understood by consider¬ 
ing a limiting case in which the emitting angles for positive 
charges are 7r/4 and Sn/A while for negative charges they are 
37r/4 and InjA. In this special limit, one can easily check 
that both jsame and jopp vanish. If the initial strength of the 
CME-like charge separation gets larger, e.g. if it is doubled 
(denoted as CESEH-2xCMEH-in-plane CSE), the magnitudes 
for both opposite-charge and same-charge correlations are be¬ 
tween those for the CME and for CESEH-CME-iin-plane CSE, 
which is shown in Pig. 4 (e). On the other hand, if we double 
the strength of the initial in-plane charge separation, e.g. in 
the case of CESEH-CMEH-2xin-plane CSE [Pig. 4 (f)], the in¬ 
plane CSE will dominate the final signal, which shows similar 
magnitudes as those for the in-plane CSE effect only [Pig. 4 
(c)]. Prom these simulations, we find that the correlation jap 
is very sensitive to CME and in-plane CSE but not sensitive to 
CESE. To test the chiral electric separation effect, we need to 
look for other observables. 


B. Two observables for CESE 

In this subsection, we propose two observables for the chi¬ 
ral electric separation effect (CESE) in Cu + Au collisions. 
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FIG. 5: (Color online) Centrality dependence of the difference be¬ 
tween opposite-charge {^opp) and same-charge {^same) correlations of 
fa/j = (cos[2(0a + for three different initial settings of AMPT 
models in Cu-l-Au collisions at yj^=200 GeV. Some points are 
slightly shifted along the horizonal axis for clearer representation. 


FIG. 6: (Color online) Centrality dependence of the difference be¬ 
tween event plane angles reconstructed from the positively charged 
hadrons and negatively charged hadrons, (I'P^ - 'P 2 I), for six initial 
settings of AMPT models in Cu-l-Au collisions at .y/j^=200 GeV. 


Let us first consider the following charge azimuthal correla¬ 
tion: 

(aH = {COS[2{(pa + (f)i3- 2'Prp)]). (2) 

We will set T'rp = 0 in the following. For a CESE + CME 
induced charge distribution [see Eig. 3 for illustration], sup¬ 
pose the positive charges fly out along the azimuthal angles 
if/c and ij/c + 7T and the negative charges fly out along -i//c and 
—ijjc + where ipc is a deformation angle due to the CESE 
effect. Thus ^same ~ (cos(4(/fc)) should be in general dif¬ 
ferent from (opp ~ 1- On the other hand, for purely dipo¬ 
lar charge distribution (due to the in-plane CSE or due to the 
CME), = (opp ~ <cos(4i/crf)), where is the dipolar an¬ 
gle {t//d = 0 for in-plane dipole and if/d - njl for out-of-plane 
dipole). Therefore, A^, the difference between ^opp and [,same, 
can be an observable for the CESE since the in-plane CSE 
or the CME does not contribute to it. In real experiments, 
the dominant background of may be the elliptic flow V 2 
because, if we turn off all the anomalous effects, ^ap ~ v\- 
Another background may be the transverse momentum con¬ 
servation [25] which causes a contribution (ap V 4 /M with 14 
the fourth harmonic flow and M the multiplicity. In these 
charge-blind V 2 and v '4 backgrounds are subtracted. However, 
there may remain other backgrounds, for example, due to the 
local charge conservation [26] (which leads to (opp v' 4 /M 
and ^same ~ 0 and thus A^ oc v '4 /M) or due to the chiral mag¬ 
netic wave induced quadrupole. Because we do not encode 
these effects into our AMPT model, the following result is 
supposed not to be related to the local charge conservation 
and chiral magnetic wave. In Eig. 5 we present the centrality 
dependence of A^ in Cu + Au collisions from three different 
initial settings in the AMPT model. One can see that the nor¬ 
mal AMPT case and the CME can not yield a visible signal for 
all centrality bins, while the result for CESEH-CME-iin-plane 
CESE case shows an increasing trend of Af from central to 


peripheral Cuh-Au collisions, where the signal becomes very 
clear in the most peripheral centrality bin. 

As discussed above, the chiral electric separation effect 
(CESE) combined with the chiral magnetic effect (CME) can 
lead to a quadrupolar charge distributions in Cu + Au col¬ 
lisions, see Eig. 3. Therefore, we propose another observ¬ 
able for CESE: AT' = <|T'+ - T'J|), where T'+ and T'^ are 
the second harmonic angles of the event planes reconstructed 
by final positively charged hadrons and negatively charged 
hadrons. If we follow the above assumption for the CESE, 
then AT* ~ 2{t(/c}. In Eig. 6 we show the centrality dependence 
of AT* for six different initial settings of the AMPT model in 
Cuh-Au collisions. As we expected, once the CESE happens, it 
leads to a sizable AT' which shows a linear dependence on the 
centrality. More importantly, for those initial settings which 
do not include the CESE, we do not observe a nonzero A*? 
for the whole centrality window. It should be mentioned that 
our results are based on the fact that in each event we know 
which side is the Au nucleus and which side is the Cu nucleus 
in the transverse plane, i.e. the direction of the electric field. 
In real experiments, it becomes very complex to identify the 
relative locations of Au and Cu and one has to reconstruct the 
first harmonic event plane in order to distinguish the Au side 
from the Cu side, which will be hopefully achieved by using 
some event plane detectors such as zero degree calorimeters 
(ZDCs) in future experiments. 


IV. SUMMARY 

In summary, we have introduced various initial charge sepa¬ 
ration effects, including the chiral magnetic effect (CME), in¬ 
plane charge separation effect (in-plane CSE), and chiral elec¬ 
tric separation effect (CESE), into the AMPT model to study 
how the final state interactions render the initial charge sepa¬ 
ration effects.To distinguish these initial effects, three possi- 
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ble observables are tested in Cu+Au collisions at ^Js^=200 
GeV. The charge azimuthal correlation (cos(0„ + <pj 3 - 2T'rp)) 
is sensitive to the CME and the in-plane CSE but not to the 
CESE. The CME results in a positive opposite-charge corre¬ 
lation and a negative same-charge correlation, while the in¬ 
plane CSE reverse their signs relative to those for the CME 
case. However, the (opp - (same is an observable that is sensi¬ 
tive to the CESE but not sensitive to CME and in-plane CSE 
which makes it be a possible observable for the detection of 
CESE. The difference between the second harmonic angles 
of event planes reconstructed by final positive-charge hadrons 
or negative-charge hadrons, (I'Tj - 'Tj|), is also sensitive to 
the CESE and is sizable in very peripheral Cuh-Au collisions. 
We thus propose that future experiments can take advantage 
of Cu + Au collisions to measure these observables in order to 
identify or distinguish different chiral and charge separation 
mechanisms. 
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